Sufficient Conditions 
for 

Fast Switching Synchronization 



Time Varying Network Topologies 



in . m 
o 
o 

<N 

D 1 Daniel J. Stilwell 1 , Erik M. Bollt 2 , D. Gray Roberson 1 

£NJ . Bradley Department of Electrical and Computer Engineering 

Virginia Polytechnic Institute and State University 
Q ' Blacksburg, VA 

u' 

2 Department of Mathematics and Computer Science 
Clarkson Unviversity 
Potsdam, NY 

February 4, 2008 

in ' 
in 

' Abstract 

| In previous work pQ, empirical evidence indicated that a time- varying network could propagate 

sufficient information to allow synchronization of the sometimes coupled oscillators, despite an in- 
stantaneously disconnected topology. We prove here that if the network of oscillators synchronizes 
for the static time-average of the topology, then the network will synchronize with the time- varying 
topology if the time-average is achieved sufficiently fast. Fast switching, fast on the time-scale of the 
coupled oscillators, overcomes the descychnronizing decoherence suggested by disconnected instan- 
taneous networks. This result agrees in spirit with that of PP where empirical evidence suggested 
that a moving averaged graph Laplacian could be used in the master-stability function analysis 
|52|. A new fast switching stability criterion here-in gives sufficiency of a fast-switching network 
leading to synchronization. Although this sufficient condition appears to be very conservative, 
it provides new insights about the requirements for synchronization when the network topology 
is time-varying. In particular, it can be shown that networks of oscillators can synchronize even 
if at every point in time the frozen-time network topology is insufficiently connected to achieve 
synchronization. 

MCS: 37A25, 37D45, 37N25, 93B05, 94C15, 05C80, 92D30 

Keywords: Nonlinear dynamics, Complex Networks, Sychronization, Controllability and observ- 
ability, self-organization, communication in complex networks. 



1 



1 Introduction 



Since Huygen's early observations of weakly coupled clock pendula 0, synchronization has been found 
in a wide variety of phenomena, ranging from biological systems that include fire flies in the forest 
Eh animal gates ^Tj, descriptions of the heart H3 [7j, and improved understanding of brain 
seizures |12j . to chemistry |15j . nonlinear optics |14| . and meteorology p2I|. See one of the many 
excellent reviews now available, including |171 1181 1191 11U1 In particular, it has been known 

for more than 20 years that chaotic oscillators can synchronize under suitable coupling |2L)1 1211 122j . 
Meanwhile, in recent years, the study of large scale, random networks has become an extremely active 
area with the advent of advances in both theory and scientific application across numerous disciplines, 
as reviewed in |231 1241 1251 1261 I27j for example. Recent investigations have sought to characterize 
how oscillator elements coupled according to a large-scale network architecture are impacted by the 
choice of architecture and corresponding spectral properties of the network, j23 EE EDI ED EH ESj- 
In particular, the master-stability function formalism [521 155j relates spectral properties of the graph 
Laplacian of the network to synchrony of supported oscillators, and this has been used in study of 
synchronization stability on arbitrary network architecture, |28| . 

Despite the very large literature to be found, the great majority of research activities have been fo- 
cused on static networks whose connectivity and coupling strengths are constant in time. For example, 
static networks are assumed for the analysis of [SU 023 US] • However, there are applications where the 
coupling strengths and even the network topology can evolve in time. Recent work such as PSIPBIETT] 
are amongst the few to consider time-dependent couplings. See also [35] in which a so-called "function 
dynamics" gives rise to networks that evolve according to a dynamical system, somewhat similarly to 
our networks. In prior work .1,, we were motivated by applications that include how a disease might 
occur in a network of agents in which the agents move, but the disease itself has its own time-scale. 
We describe a competition of two time-scales. Said plainly, the disease has a natural typical incuba- 
tion rate, and a natural infections rate (for example in a Susceptible— > Exposed— ^Infected— > Recovered, 
SEIR, model) , so if a susceptible agent does not come in contact with an infected agent in the disease 
time scale, then there should be no new infection. Mathematically, we constructed a "moving neigh- 
borhood network" (MNN), a network of agents which move ergodically, and connect when in close 
proximity to each other. Such a network was shown to lead to a Poisson distributed degree distribu- 
tion, instantaneously, and hence the network was typically instantaneously disconnected. It consists 
of typically many small subcomponents at each instant. Such a description alone would suggest that 
there could be no global synchronization of the oscillators carried by each agent which are coupled ac- 
cording to the disconnnected network. However, it was found that if the agents move quickly enough, 
then roughly described, in any recent time window, a given agent might be likely to have had some 
amount of coupling to most other agents. It turns out that for fast enough moving agents, these ran- 
dom time varying connections were enough to overcome even chaotic oscillators' sensitive dependence 
tendency to drift apart asynchronysly. We formalized this idea by introducing a new description of 
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the connectivity, a "moving averaged" graph Laplacian. We showed empirically that the spectrum of 
this construction works quite well together with the master stability formalism to accurately predict 
synchronization. 

Besides our original motivation in mathematical epidemiology, it can be argued that this work has 
strong connections to ad hoc communication systems and control systems on time varying networks. 
Fundamental connections between chaotic oscillations and proof of synchronization through symbolic 
dynamics |391 14U| and control |41l 1421 I43| and even definition of chaos through symbolic dynamics 
suggest this work is rooted in a description of information flow in the network. 

Coordinated control for platoons of autonomous vehicles can also be addressed using network 
concepts |58l I6U1 159*] . Each vehicle is represented by a node, and communication or mutual sensing 
is represented by connections between nodes. In the average position of a platoon of vehicles is 
regulated, and the graph Laplacian is used to describe communication between vehicles. It is shown 
that the spectrum of the graph Laplacian can be used to indicate stability of the controlled system. 
As pointed out in [HO], the use of a graph Laplacian is not entirely common since it appears naturally 
for only a limited class of control objectives. 

These considerations have led us in this work to consider a simplified version of the moving agents 
of our MNN model. Considering certain time varying coupled network architectures, we can now make 
rigorous but sufficient statements concerning fast switching, and we use mathematical machinery not 
so far typically used in the synchronization community. The main result of this work comes from the 
fields of switched systems, and specifically builds on the concept of fast switching. Switched systems 
are a class of systems whose coefficients undergo abrupt change, for example, consider the linear state 
equation 



where pit) : M. i— > Z + is a switching sequence that selects elements from a family of matrix-valued 
coefficients = {A%, A2, ■ ■ •}■ When each element of is Hurwitz, stability of Q is guaranteed if 
p{t) switches sufficiently slowly. Further restrictions on elements of O, such as existence of a common 
Lyapunov function, can guarantee stability for arbitrary switching functions, including those that are 
not slow. An excellent overview of the field of switched systems and control is presented in |44j and 
in the book |45j . 

Even when the elements of are not all Hurwitz, stability of Q is still possible, although the 
class of switching functions is further restricted. For example, in jIBf a stabilizing switching sequence 
is determined by selecting elements of based on the location of the state x(t) in the state space. 
This is essentially a form of state feedback. 

When no elements of are Hurwitz, which is the case that is considered herein, stability of ([T]) 
can sometimes be guaranteed if the switching sequence is sufficiently fast. Loosely speaking, it can be 
shown that 



x(t) = A p(t) x(t) 




Xit) = A p{t/e)x{t) 
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is asymptotically stable if the time-average 

1 f t+T 

- J A p{u) da 

is Hurwitz for all t, and if e is sufficiently small. This fact has been established in (471 1481 149j for 
several classes of linear systems related to (J2J. Similar results have been presented in |5U1 151j for 
classes of nonautonomous nonlinear systems where time is parameterized by t/e as in In this 
case, stability of a specific average system implies stability of the original system if e is sufficiently 
small. In addition, this work requires the existence of a Lyapunov function that is related to a certain 
average of the system but which is not a function of time. This requirement is too restrictive for 
the class of linear time-varying systems considered herein. A new fast switching stability condition, 
presented in Section |31 is derived in order to assess local stability of networked oscillators about the 
synchronization manifold. 



2 Preliminaries 

We consider a network of coupled oscillators consisting of r identical oscillators, 

r 

± i (t)=f(x l (t))+aBY,hj(t)x j (t), i = l,...,r (3) 
i=i 

where X{(t) € M n is the state of oscillator i, B E M nxn , and the scalar a is a control variable that sets 
the coupling strength between oscillators. The scalars lij{t) are elements of the graph Laplacian of 
the network graph and describe the interconnections between individual oscillators. Let G(t) be the 
time- varying graph consisting of r vertices Vi together with a set of ordered pairs of vertices {vi , Vj } 
that define the edges of the graph. In this work, we assume that {vi,Vi} € G(t) for i = 1, ... ,r. Let 
G(t) be the r x r adjacency matrix corresponding to G(t), then Gij(t) = 1 if {vi, Vj} is an edge of the 
graph at time t and Gij(t) = otherwise. The graph Laplacian is defined as 

L(t) = diag(d(t)) - G(t) (4) 

where the i th element of d(t) € W is the number of vertices that vertex % is connected to, including 
itself. 

Synchronization can be assessed by examining local asymptotic stability of the oscillators along the 
synchronization manifold. Linearizing each oscillator (J2J) about the trajectory x°(t), which is assumed 
to be on the synchronization manifold, yields 

r 

Zi(t) = F(t)zt(t) + aBj2 l ij(t)zj(t) (5) 

3=1 

where Z{(t) = Xi(t) — x°(t), and F(t) = Df evaluated at x°(t). Let L(t) be the r x r matrix composed 
of entries hj(t), then the system of linearized coupled oscillators is written 

z(t) = (J r ® F(t) + a(I n ®B)(L® I r )) z(t) 

(6) 

= (J r ® F(t) +aL®B) z(t) 
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where '(8)' is the Kronecker product and z(t) = [zf(t), . . . , zJ(t)] T . Standard properties of the Kro- 
necker product are utilized here and in the sequel, including: for conformable matrices A, B, C, and 
D, (A (g) B){C ®D) = AC ® BD. Notation throughout is standard, and we assume that || • || refers to 
an induced norm. 

It has been shown in |52l 15*5] that the linearized set of oscillators (JHJ) can be decomposed into two 
components: one that evolves along the synchronization manifold and another that evolves transverse 
to the synchronization manifold. If the latter component is asymptotically stable, then the set of 
oscillators will synchronize. 

The claimed decomposition is achieved using a Schur transformation. We briefly describe the 
decomposition since it plays a central role in our assessment of synchronization under time-varying 
network connections. Let P € K nxn be a unitary matrix such that U = P~ 1 LP where U is upper 
traigular. The eigenvalues Ai,...,A r of L appear on the main diagonal of U. Then a change of 
variables £(t) = (P ® I n )~ z(t) yields 

((t) = (P ® In)' 1 (J r ® F(t) +oL®B){P® I n ) f (i) 

= (I r ® F(t) + aP~ l LP ® B) e(t) (7) 
= (J r ® F(t) + all B) 

Due to the block-diagonal structure of I r <g> F(t) and the upper triangular structure of U, stability of 
(JZJ) is equivalent to stability of the subsystems, 

£ i (t) = (F(t)+a\ i B)£ i (t), i = l,...,r (8) 

where Ai, . . . , X r are the eigenvalues of L. Note that since the row sums of L are zero, the spectrum 
of L contains at least one zero eigenvalue, which we assign Ai = 0. Thus 

6(t) = F(t)ei(o 

evolves on the along the synchronization manifold, while (jSJ) with i = 2, . . . , r evolves transverse to 
the synchronization manifold [52]. Since the oscillators are assumed identical, the (identity) synchro- 
nization manifold is invariant for all couplings, the question being its stability. The set of coupled 
oscillators will synchronize if the synchronization manifold is stable, if © with % = 2, . . . , r is asymp- 
totically stable. 

3 Main Result 

For a given static network, the master stability function characterizes values of a for which a set of 
coupled oscillators (jBJ) synchronizes |52l 1531 154] . The graph Laplacian matrix L has n eigenvalues, 
which we label, 

= Aq < • • • < A n _i = Xmax- (9) 
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The stability question reduces by linear perturbation analysis to a constraint upon the eigenvalues of 
the Laplacian: 

crXi G (0,1,0.2) Vi = 1, . . . , n — 1, (10) 

where ct\, 02 are given by the master stability function (MSF), a property of the oscillator equations. 
For a small, synchronization is unstable if a \\ < a±; as a is increased, instability arises when, 



O-Xmax > OL 2 . (11) 



By algebraic manipulation of (fTU)l . if, 



<?*=:0, (12) 
Ai a\ 

then there is a coupling parameter, a s , that will stabilize the synchronized state. For some networks, no 
value of a satisfies (|10|> . In particular, since the multiplicity of the zero eigenvalue defines the number 
of completely reducible subcomponents, if Ai = 0, the network is not connected, and synchronization is 
not stable. However, even when Ai > 0, if the spread of eigenvalues is too great, then synchronization 
may still not be achievable. 

For the case of a time- varying network topology, represented by Lit), our principal contribution 
is to show that the network can synchronize even if the static network L = Lit) for any frozen value 
of t is insufficient to support synchronization. Specifically, we show that the time-average of L(t), not 
the frozen values of Lit), is an indicator of synchronization. If the time-average of L(t) is sufficient 
to support synchronization, then the time-varying network will synchronize if the time-average is 
achieved sufficiently fast. 

Theorem 3.1 Suppose a set of coupled oscillators with linearized dynamics 

z s (t) = (l r ®F(t) + aL®B)z s {t) (13) 

achieves synchronization. Then there exists a positive scalar e* such that the set of oscillators with 
linearized dynamics 

z a (t) = (I r ® F(t) + aL(t/e) ® B) z a {t) (14) 

and time-varying network connections L(t) achieves synchronization for all fixed < e < e* if L(t) 
satisfies 



1 

f 

and the column sums of L(t) are all zero for all t. 



t+T 

L(a)da = L (15) 



Remark 3.2 Since L(t) represents a time-varying network, we may assume that for each value oft, 
L{t) is a graph Laplacian as defined in Q. Thus the time-average L in (|15|) is not a graph Laplacian. 
In other words, L does not necessarily correspond to a particular network topology and arises only as 
the time-average of L(t). However, L does inherit the zero row and column sum property of L(t). 
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A preliminary lemma is required to prove Theorem 13.11 the proof of which appears in the Appendix. 



Lemma 3.3 Suppose the matrix-valued function E(t) is such that 



1 

T 



t+T 



E(a)da = E 



for all t and 



(16) 
(17) 



x(t) = (A(t) + E)x(t), x(t Q ) =x Q , t> t 
is uniformly exponentially stable. Then there exists e* > such that for all fixed e € (0, £*), 

z(t) = (A(t) + E(t/e))z(t), z{t ) =z a , t> t Q (18) 

is uniformly exponentially stable. 

Proof of Theorem \S.l\ 

First we show that the Schur transformation that decomposes the set of oscillators (|13|) with static L 
also induces a similar decomposition for (|14[) with time- varying L(t). Then we apply Lemma 13.31 to 
show that the modes of the system that evolve transverse to the synchronization manifold are stable 
if e is sufficiently small. 

Let P £ R rxr be a unitary matrix such that U = P~ 1 LP where 



U 



Ui 
0( r -l)xl U 2 _ 



r i)x(r i) j g U pp er triangular. Without loss of generality, 

1/r, y / l/r 



is the Schur transformation of L, and U2 € 

we have assumed that the left-most column of P is the unity norm eigenvector 
corresponding to a zero eigenvalue. The change of variables £ s (t) = (P ® I)~ 1 z s (t) yields the decom 
position £ s = £s2] T where £ s i G M. n , £ S 2 € M™^ -1 ), and ^ S 2 satisfies 



U(t) = (J r _i ® F(t) + 0-U2 ® B) U(t) 

As discussed in Sectional (|19|) is asymptotically stable by hypothesis. 

We now consider the same change of variables applied to (|14jl. First, note that 



(19) 



u(t) = p^nfyp 



Ut(t) 

0( r -l)xl U^it) 



since the column sums for L(t) are zero for all t. The change of variables £ a (i) = (P <8> I)~ 1 z a (t) 
yields the decomposition ^ a = [£ a i>£a2] T where € W 1 evolves along the synchronization manifold 
and £ a 2 G R n ( r_1 ) evolves transverse to the synchronization manifold. To verify that the oscillators 
synchronize, it is sufficient to show that 



= (Ir-i ® F(t) + aU 2 (t/e) (8> B) £ a2 (t) 



(20) 
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is asymptotically stable when e is sufficiently small. Since 

U = P~ l LP 

r-t+T 



1 

T 
1 

f 



P^LMPda 



t+T 

U{a)da 



we conclude that ^ ^ 

U 2 = ^j t + U 2 (a)da (21) 

Thus the desired result is obtained by direct application of Lemma 1331 along with (|19[). (|20|). and (|21|). 
■ 

4 Illustration 

To illustrate fast switching concepts applied to synchronization of a set of oscillators, we consider a 
set of r Rossler attractors 

r 

±i(t) = -Viit) - Zi(t) - <jy2 l ij(t/e)xj(t) 

(22) 

jfi(t) = Xi(t) + ayi(t) 
zj(t) = b + Zi(t)(xi(t) - c) 

where i = 1, . . . ,r, a = 0.165, 5 = 0.2, c = 10, and er = 0.3. Oscillators are coupled through the Xi 
variables via kj(t). Coupling between subsystems (nodes) is defined by a time- varying graph G(t), 
with corresponding adjacency matrix G(t). The graph Laplacian L(t), with entries hj(t) is defined as 
in 

For the purposes of illustration, we choose a set of five graphs and corresponding adjacency matrices 
G±, . . . , Gr5 with the property that none of them are fully connected. That is, each graph contains 
pairs of nodes that do not have a path between them. However, the union of vertices over all five 
graphs yields a fully connected graph with the longest path between nodes containing no more that 
two other nodes. All five graphs and the union of graph vertices are shown in Figure ^ 

A simple strategy is chosen for switching among graph Laplacians associated with the set of graphs. 
We choose the T-periodic L(t) defined over one period by 

5 

L{t) = ^ LiX[(j-i)T/5, iT/5) (t) 

8=1 
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where X[u,t2)(t) is the indicator function with support [t±, t 2 )- The time-average of L(t) is 

i r T 

L = 7TJ L{t/£)dt 

1 5° (23) 

= |S> 

i=l 

Toward computing the upper bound for e given by ([35 j). the set of coupled oscillators ((2*2*|) with 
coupling defined by (|23|) are integrated. The cc-coordinate for each oscillator is shown in Figure 
12 The oscillators clearly synchronize. Asymptotic stability of the oscillators with respect to the 
synchronization manifold is suggested by plotting the sum-square deviation of the states 

r 

£(xi(t) - p x {t)f + ( yi (t) - p y (t)f + (zi(t) - p z {t)f (24) 
i=l 

about the averages 

1 r 

i=i 

where /U^(t) and /^(i) are defined similarly. Approximately exponential decay of (|24|) is evident in 
Figure 01 

The linear time-varying system © corresponding to the set of coupled Rossler attractors is com- 
puted from the Jacobian of the right-hand side of 1)221) evaluated at the solutions shown in Figure 

121 

As described in the proof of Lemma l.3.31 a Schur transformation U that diagonalizes L is computed 
and used as a state transformation to decompose the linear time- varying system (jf)J) into a component 
that evolves along the synchronization manifold and another component that evolves transverse to 
the synchronization manifold. The upper bound for e given in Theorem 13.11 is computed from the 
component of the linear system that evolves transverse to the synchronization manifold, 

la(t) = (ir-i ® F(t) + aU 2 ® B)Cas(t) 

We now estimate the constants a, p, rj, and \x needed to compute the right-hand side of (|33)) in the 
proof of Lemma 13.31 (see Appendix). This is used to compute an maximum value of e. The constant 
a is computed from ()28|) . while the transition matrix is computed from 

$(t,r) = (J r _i®F(t) +aU 2 ®B)<f>(t,T), $(r,r) = / 

The norm of the transition matrix ||<&(t, t)|| is shown in Figure HJ The initial time r is chosen to be 
40 seconds to ensure that the states of (1221) are reasonably close the the synchronization manifold. An 
upper bound that satisfies r)|| < ^e~ x<yt ~ T ' > is also shown in Figure 0J The coefficients p, p and 
r] in ()32|) are computed from 7 and A when evaluating the right-hand side of (|35|) . Choosing T = 1, 
the right-hand side of Q35|) is evaluated for this example, and we determine that the set of couple 
oscillators will synchronize if e < 3.3 x 10~ 7 . This shows that our bound is exceedingly conservative. 
For example, empirically the oscillators will synchronize with e = 1, as shown in Figure |SJ 
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(d) (e) (f) 

Figure 1: (a)-(e) are graphs G\ through G5, respectively, while (f) is the union of graphs. 
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Figure 3: Root-mean-squared error in (j2D) for the set of coupled Rossler attractors using the average 
network L. 
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Figure 4: Norm of the transition matrix <J>(i, r) along with an exponentially decaying upper bound. 
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Figure 5: The x-coordinate for the set of coupled Rossler attractors using the switched network where 

8 = 1. 
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6 Appendix 



Proof of Lemma 13, 31 

Since (|17|) is uniformly exponentially stable, there exists a symmetric matrix function Q(t) and 
positive scalars rj, p, and p, such that the Lyapunov function 

v(x(t),t) = x T (t)Q(t)x(t) 

satisfies 

r]\\x(t)\\ 2 <v(x(t),t) < p\\x(t)\\ 2 (25) 

j t v{x{t),t) < -p\\x{t)\\ 2 (26) 

for all t. To establish uniform exponential stability of we will show that v(z(t),t) is also a 

Lyapunov function for (|18|) if e is sufficiently small. This claim is achieved by showing that for 
sufficiently small values of e, 

Av(z, t + eT, t) = v(z(t + eT),t + eT) - v(z(t),t) (27) 

is negative definite for all t. Expanding (|27jl yields 

Av(z, t + eT, t) = z T (t + eT)Q(t + eT)z(t + eT) - z T (t)Q(t)z(t) 

= z T (t) U T E {t + eT, t)Q(eT + t)$ E (t + eT, t) - Q(t)) z(t) 



where & E (t,r) is the transition matrix corresponding to A{t) + E(t/e). Similarly denoting the tran- 
sition matrix for A{t) + E as § E {t,r), we define 

H{t + eT,t) = <5> E {t + eT,t) - <5> E {t + eT,t) 

r-t+eT 



/t+el 
A{oi) + E{a/e)d<j 

00 rt+eT ra x ra i - 1 

+ V / A(ai) + E(a 1 /e) •••/ A(<7j) + E^fe)^ ■ ■ ■ da t 

■r, Jt Jt Jt 



It Jt Jt 

rt+eT 00 rt+eT ra 1 ra i _ 1 

-I- A(a 1 ) + Edcr -V / A{a x ) + E / •••/ A(&i) + Edoi • • • do\ 

Jt i=2 Jt Jt Jt 

By hypothesis, 



rt+eT 

j E(a/e)da = eTE 



which implies that 

x rt+eT ra x rvi-\ 



rt+ei ra\ r^i-i 

H(t + eT,t) = J2 A (°i) + Eiax/e) / • • • / Afa) + E(ai/e)dai ■ ■ ■ da x 

i=2 Jt Jt Jt 

00 rt+eT ra± rcti-\ 

- V / A(ai) + E ■■■ A{ai) + Edoi ■■■da 1 

■ o Jt Jt Jt 
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Defining 



a = sup(max(\\A(t)+E\\,\\A(t) + E(t/e)\\)) (28) 
t>o 



a bound for H(t + eT, t) is computed 

\\H(t + eT, t)\\ < 2{e eTa - 1 - eTaj (29) 
Noting that <&e = &e + H, Av is expressed 

Av(z, t + eT, t) = z T (t) ($|(t + eT, t)Q(t + eT)$ E (t + eT, t) - Q(tfj z(t) 

+ z T {t) ($|(t + eT , i)Q(* + eT)fl"(t + eT, t) + tf T (t + eT, t)Q(t + eT)^(t + eT, t) 
+ H T (t + eT, t)Q(t + eT, t)H(t + eT, t)\ z(t) 

(30) 

The task now is to compute an upper bound for Av(z, t + eT, t) and show that this bound is negative 
if e is sufficiently small. Several well-known relationships that are consequences of (|25|). (|26j) . and 
uniform exponential stability of (|17j) are utilized (see for example JSH] pages 101 and 117, or |S7j, page 
202). Namely, 

\\Q(t)\\<P (31) 
|[^(Mo)||<y^e-£ (t - to) (32) 
v(x{t),t) < e- E p {t - to) v{x{t ), t ) (33) 

for t >t Q . 

To compute an upper bound for the first term on the right-hand side of (|30[) we note that if 
x(t) = z{t) is chosen as the initial condition of (|T7|) at time t, then 



z 1 (t) ($|(i + eT, t)Q(t + eT)$ E (t + eT, t) - Q(t)j z(t) = v(x(t + eT),t + eT) - v(x(t),t) 
Prom (jS3) and (HSJ), 

v{x(t + eT),t + eT) - v(x(t),t) < ( e ~^ T/p - l)v(x{t),t) 

< p ( e -^ T /p -l)\\ x (t)f 

Thus, 

z T {t) (<J?|(i + eT, t)Q(t + eT)$ E (t + eT, t) - Q(t)) z(t) < p(e~^ T /P - l)\\z(t)\\ 2 (34) 
Combining (|29|). (|3T]). (|3*2*)) . and (|3*1|) yields the desired upper bound 



Au(z,t + eT,t) < (p(e-^ T/p - 1) + 4 / o( v / pA?e~^ r )(e £Ta - 1 - eTa) + 4p(e eTa - 1 - eTa)- ! | : 



(35) 
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Defining the continuously differentiable function g{e,x) to be the right-hand side of (|35[). it can be 
shown that g(0, z) = and -§^g(0, z) = — fj,T\\z\\ 2 < 0. Thus since g(e, z) — > oo as e — > co, there exists 
e* such that g(e*,z) = and z) < for all e G (0, £*) and z^O. Thus Aw(z, t + eT,t) < for all 
£ G (0,e*) and z / 0. 

To show that negative-definiteness of Av(z,t + eT, t) is sufficient to establish stability of (j!8[) . Choose 
e and 7 > that satisfy 

Av(z,t + eT,t ) = v(z(t + eT),t + eT) -v(z(t ),t ) < - 7 ||z(t )|| 2 

for all t Q . From (|25|). w(z(t ),t ) < /9||z(t ) |[ 2 , which implies that 

+ eT),t + eT) - v(z(t ),t ) < -(>y/p)v(z(t ), t Q ) 

Thus 

v{z{t Q + eT),t + eT) < (1 - i/p)v{z{t ),t ) 
Repeating this argument yields 

v{z(t + keT),t + keT) < (1 - ^ / P ) k v{z{t ), t Q ) 

for any positive integer k. Thus v(z(t +keT),t +keT) — > as — > oo which implies that z(t Q +keT) — > 
as A; — > oo. Since the limiting behavior is valid for any t , uniform exponential stability of is 
established. ■ 
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